.ry 
r.  3.  Naval  Postgraduate  School 
Monterey,  California 


SOME  EXAMPLES  OF 
EXPONENTIAL  REFRACTIVITY  DISTRIBUTIONS 


by 

Charles  G.  Waldron 
// 

Lieutenant,  United  States  Naval  Reserve 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 

METEOROLOGY 

United  States  Naval  Postgraduate  School 
Monterey,"  California 

19  6  0 


^>y&°3 


School  . 


SOME  EXAMPLES  OF 
EXPONENTIAL  REFRACTIVITY  DISTRIBUTIONS 
by 
Charles  G.  Waldron 


This  work  is  accepted  as  fulfilling 

the  thesis  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 

METEOROLOGY 

from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 

Recently,  exponential  models  of  refractivity  have  been 
devised  by  several  investigators  to  replace  the  linear,  "standard" 
refractive  atmosphere.  In  this  study,  average  values  of  N  are 
evaluated  from  the  data  at  the  mandatory  radiosonde  levels  up 
to  25  mb,  and  mean  coefficients  c^  satisfying  A/U/=/V.  a"c^2cM~z;J 
have  been  evaluated  for  a  group  of  eight  mid-latitude  stations, 
and  for  the  entire  winter  season  of  1958-59.   The  coefficients 
c^,  plotted  as  functions  of  the  height  of  the  i-th  layer- 
center,  possess  certain  characteristic  features  in  common  at  all 
stations. 

The  writer  wishes  to  express  his  appreciation  for  the 
assistance  and  enoouragement  given  him  by  Professor  Frank  L. 
Martin  of  the  U.  S.  Naval  Postgraduate  School  in  this  investiga- 
tion. 
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1.   Introduction 

The  refractive  index  n  is  defined  as  the  ratio  of  the 
speed  of  an  electromagnetic  wave  in  a  vacuum  to  the  wave  speed 
in  the  medium  concerned.   However,  this  particular  form  of  the 
refractive  index  is  not  ve  of 

ray-paths;  hence  the  modified  refractive  index,  N  is  normally 
used. 

The  refractivity  N,  which  is  equal  to  (n-1)  x  10  ,  may 
be  obtained  from  radiosonde  observations  using  the  equation 


N    ,    5.   [p  i  -£_£)   =  '  (-n-i  )  x  jo 


6 

(1) 


where  T  is  the  temperature  in  degrees  Kelvin,  p  is  the 

pressure  in  millibars,  e  is  the  vapor  pressure  in  mbs,  and 

A  and  &  are  constants. 

Using  the  Smith  &  Weintraub  [4]  values  of  these 

constants,  Equation  (1)  may  be  written: 

a/     1  1    Q     P  373Z5G-K- 

A/     =     n^  +  — — ^ (2) 

A  "standard"  distribution  of  N  in  the  vertical  may  be 
obtained  by  averaging  over  a  sufficient  number  of  individual 
soundings.   Bean  and  Thayer  [lj  have  suggested  that  its 
curve  can  be  fitted  by  an  exponential  formula.   The  model 
which  has  recently  been  adopted  by  the  National  Bureau  of 
Standards  is  of  form: 

N    ~-     A/s   e  ~e  Z  P) 

where  N  is  the  refractive  index  at  any  altitude  , 


Ng  is  the  mean  surface  value  of  the  refract] 

z  is  the  height  above  the  surface, 

c  is  an  empirical  constant,  affording  "best  fit"  in 
Equation  (3). 

Equation  (3)  is  used  only  with  the  mean  values  of  N 

averaged  over  a  long  enough  period  to  eliminate  all  transitory 

trapping  layers.   The  study  by  Bean  and  Thayer  flj  also  suggested 

further  study  by  obtaining  individual  values  of  Ns  and  c  for 

varying  climatic  regions.   This  was  the  objective  of  the 

first  approach  of  this  author  to  the  problem.   Later,  however, 

it  was  decided  to  make  a  finer-scale  determination  of  the 

parameter  c  by  determining  a  mean-value  c^  for  each  layer 

between  mandatory  levels.   Thus  a  working  equation  of  form 

^U>    '  "\    e~C;C*'2i)   ,   2,  <  Z   <  ZltJ  (4) 

was  employed.   Here  c^  is  temporarily  assumed  to  be  constant 
in  the  i-th  layer,  which  extends  from  z^  to  z^+^  and  Nj[  is 
the  value  of  N  at  z  =  z^     Note  finally  that  c£  may  be 
obtained  by  differentiation  of  (4)  as 

C  ■     --     -   -i-  $-£±  t^ 


2.   Data  Processing 

The  stations  utilized  in  this  study  were  selected  to 
represent  a  wide  variety  of  climatic  regions.   The  cj- values 
were  determined  from  the  three-month  means  of  pressure,  tempera- 
ture, and  moisture  content,  for  all  the  mandatory  levels  at 
these  stations  during  the  winter  of  1958-59.   These  data  were 
obtained  from  the  radiosonde  checked  data  found  in  the 
Northern  Hemisphere  Data  Tabulations  Daily  Bulletins  [2J  . 

The  use  of  data  obtained  from  radiosondes  introduced  some 
uncertainty  of  the  relative  humidity  when  "motorboating"  was 
reported.   However,  motorboating  usually  occurred  near  400  mb , 
and  there  the  error  introduced  by  omitting  "motorboating" 
reports  from  the  moisture-  averaging  process  was  found  to  be 
^0.1  mb.   The  twice—daily  observations  taken  at  0000  and 
1200  Greenwich  times  were  averaged.   This  served  to  filter 
out  the  diurnal  variation,  and  to  make  the  stations,  some  of 
which  were  separated  by  as  much  as  six  hours  time-difference, 
more  compatible.   Two  methods  were  used  in  obtaining  the 
values  utilized  for  the  determination  of  the  surface  refractive 
index  data.   When  the  surface  report  was  part  of  the  published 
radiosonde,  this  was  the  preferred  data-source,   V7hen  this 
radiosonde  level  was  not  available,  it  was  necessary  to  employ 
the  surface  observation  for  the  station.   However  the  surface 
observation  was  recorded  only  once  each  day,  at  1200Z,   and 
this  would  result  in  a  comparison  of  a  once-daily  observation 


against  a  twice  daily  one  when  going  from  the  surface  to  850 
millibars.   As  a  test,  the  surface  three-month  mean  N-value  for 
Washington  was  computed  using  both  methods  and  a  difference  of 
only  one  N-unit  in  300  was  noted.   Since  Washington  had  one  of 
the  largest  surface  diurnal  variations  ." , 

this  indicated  that  the  once-a-day  surface  observation  reports 
could  be  substituted  for  the  two-a-day  radiosonde  surface 
reports  without  serious  error. 

For  the  actual  determination  of  mean  values,  the  standard 
levels  from  the  surface  up  to  25  mbs  were  used.   The  twice-daily 
observations  at  these  levels  were  totaled  and  then  divided  by 
the  total  number  of  entries.   All  the  calculations  were  done  on 
either  tape  adding  machines  or  hand  calculators.   A  maximum  of 
180  observations  were  used  for  each  lev^l,  pn 

.  The  mean  vapor  pressure  was  obtained  by  determining 
the  vapor  pressure  from  each  observation  using  the  temperature 
and  relative  humidity  data,  and  finding  the  mean  of  these.   For 
the  majority  of  the  stations,  moisture  was  recorded  only  up  to 
the  400-mb  level,  at  which  level  the  vapor  pressure  was  found  to 
be  less  than  1  mb.   The  only  station  that  had  moisture  data  for 
the  300-mb  level  showed  the  mean  vapor  pressure  to  be  less  than 
one-tenth  of  a  mb.   This  showed  that  very  little  accuracy  would 
be  lost  by  not  including  moisture  considerations  above  400  mbs. 

After  the  three-month  mean  values  of  pressure,  temperature, 
and  vapor  pressure  were  obtained  for  the  mandatory  levels  of 
each  station,  the  three-month  mean  Nj[  values  were  determined 


for  these  levels  using  Equation  (2).   These  values  and  the 
appropriate  pressure  levels  are  given  in  Table  1. 

The  C£  -values  for  each  layer  were  obtained  using  the 
refractivity  equation  in  the  form: 

*~  #" 

C>   =  -  2.  ,-zl  (6) 

where  z^+^  -z^  is  the  average-value  of  the  thickness  in  the 

standard  layers  employed  in  this  study. 

The  values  of  c^  for  each  station  were  then  plotted  on  two 
cycle  semi-logarithmic  paper  versus  the  midpoint-height  of  each 
layer  as  ordinate.   These  midpoints  were  taken  from  the  U.  S. 
Standard  Atmosphere  elevations  for  the  stated  pressure.  A 
discussion  of  these  results  and  of  some  applications  are  given 
in  the  following  sections. 
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3.   Discussion  and  Results 
(a)  Typical  Model 

By  Equation  (5),  Cj[  is  a  good  indication  of  the  vertical 
distribution  of  N.   An  increasing  value  of  ci  indicates  a 
larger  rate  of  decrease  of  N  with  height.   Most  stations 
studied  showed  two  main  values  of  c^  in  the  vertical.   A 
nearly  constant  value  was  indicated  for  the  troposphere,  and  a 
somewhat  larger  value  shown  for  the  stratosphere.   These 
values  were  separated  at  most  stations  by  a  characteristic 
feature,  which  will  be  referred  to  as  a  dual  discontinuity. 
The  dual  discontinuity  is  evidently  associated  with  the  tropo- 
pause;  it  can  be  described  as  having  a  minimum  value  of  c^ 
and  immediately  above  it  a  pronounced  increase  to  a  maximum 
value.   A  model  of  c^  vs  height  which  depicts  schematically 
the  association  between  the  dual  discontinuity  and  the 
related  temperature  sounding  is  shown  in  Figure  1,  which  may 
be  regarded  as  typical  of  most  mid-latitude  stations. 

In  order  to  understand  the  association  depicted  in 
Figure  1,  consider  N  as  given  by  Equation  (2)  with  e=0.   Then 

N   -    77c(f  )  (7) 

Partial  differentiation  of  N  in  (7)  with  respect  to  z, 
and  use  of  the  hydrostatic  equation  leads  to 
d  N     _  _  7  7.6  P  (  v\ 

where   K    is  the  autoconvective  lapse  rate,  and   #   is  the 
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Figure  1,   Schematic  relationship  between  the-  radiosonde 
and  the  associated  refract ivity  profile,  for  mid- 
latitude  stations. 
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atmospheric  lapse  rate.   By  use  of  Equation  (5),  it  follows 
that 

C;  =  ^r   (  fA-  0  (8) 

Hence  c^  is  relatively  large  in  a  thermally  stable  layer 
and  relatively  small  in  a  thermally  unstable  layer.   Thus  if 
we  disregard  moisture,  any  minimum  value  of  c^  is  associated 
with  a  relatively  steep  lapse  in  the  layer,  and  a  maximum 
value  with  a  relatively  stable  lapse.   The  discussion  of  the 
individual  stations,  which  follows,  will  point  out  some 
deviations  from  the  typical  model,  Figure  1. 

One  can  also  qualitatively  estimate  what  the  effect  of 
moisture  lapse  would  be  on  c^.   Empirically  one  knows  that  a 
large  negative  value  of  dhJ/Az     may  be  associated  with  a 
strong  hydrolapse  and  this  would  cause  c^  to  be  large 
[Equation  (5)J   .   On  the  other  hand,  if  the  temperature 
lapse  is  not  abnormal,  c^  may  be  a  minimum  as  a  result  of  a 
small  lapse  (or  possibly  an  increase  with  height)  of  water 
vapor.   Moisture  effects  may  be  associated  with  some  of  the 
deviations  of  the  individual  c^  graphs  from  the  ideal  one 
of  Figure  1. 
(b)  Adak 

The  tropospheric  distribution  of  c^  vs  height, 
Figure  2,  showed  a  gradual  decrease  of  c^  with  increasing 
height  up  to  9  km.   Above  this  level,  a  pronounced  dual  dis- 
continuity was  found,  indicating  the  presence  of  a  tropopause 
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layer.   A  second  dual  discontinuity  of  lesser  magnitude  was 
noted  above  the  first.   This  was  taken  as  an  indication  of  a 
multiple  tropopause,  which  is  in  keeping  with  the  model  of 
Figure  1. 

(c)  Denver 

The  graph  for  Denver,  Figure  3,  shows  a  smaller  value  of 
C£  in  the  troposphere  than  in  the  stratosphere.   There  is  a 
sharp  increase  in  Cj_  to  a  maximum  value  near  13  km,  with  no 
definite  dual-discontinuity  minimum  below. 

(d)  Ship  "P",  Tatoosh  Island,  Kodiak,  Washington,  D.C. 
The  graphical  results  for  these  stations  are  shown  in 

Figures  4,  5,  6,  7,  respectively.   These  graphs  have  been 
grouped  together  because  they  all  have  features  closely  re- 
sembling the  idealized  model  of  Figure  1.   The  tropospheric 
cj_  decreases  slightly  with  elevation,  and  is  followed  by  a 
faster  rate  of  decrease  ending  in  the  region  of  8"  10  km. 
The  stratospheric  c^  then  conforms  rather  closely  to  the 
schematic  model.   Kodiak  shows  a  lesser  tendency  to  have  the 
pronounced  minimum  c^  as  part  of  its  dual  discontinuity. 
This  might  well  be  a  characteristic  of  its  higher  latitude, 
together  with  a  greater  variability  of  tropopause  height. 

(e)  Ship  "N" 

The  graphical  results  for  this  station,  Figure  8,  shows 
a  greater  decrease  of  tropospheric  c^  with  elevation  than  is 
shown  in  the  idealized  model.   This  is  due  primarily  to  the 
greater  moisture  content,  and  consequently  greater  N-value,  in 
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low  leveLs  at  this  subtropical  maritime  station.   Thus  the  low- 
level  values  oi:  e^  are  larger  than  those  for  most  other  stations, 

although  the  minimum  just  below  the  tropopause  has  nearly  the 

-5 
same  value  (about  11.5  |c  10   per  meter). 

(£)   Lake  Charles 

The  graph  for  (CSke  Charles,  Figure  9,  shows  the  tropo- 

spheric  c^  -value  decreasing  to  a  minimum  which  extends  from 

6  km  to  1.0  km.  Above  this  layer,  c^  increases    a     l" 

constant  maximum  value  in  the  stratosphere.   Tht  "ant 


cj        range  6-10  km  is  a  characteristic  of  variable  tropo- 
t  Lght,  during  the  observation  period  (c.f.  Kodiak, 

F2  ,jre  6). 

- 

(g)   Coi.'.j;  ;site  of  all  stations 

The  composite  graph  of  all  stations,  Figure  10,  conforms 
closely  to  the  idealized  model.   The  ci  -value  decreases 
slowly  with  height  in  the  troposphere,  and  is  separated  from 
the  larger  stratospheric  cj_    -values  by  the  dual  discontinuity. 
A  secondary  tropopause  is  also  indicated  by  the  less  pronounced 
upper  dual  discontinuity  in  the  neighborhood  of  15  km, 
within  the  stratosphere. 
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4,   Othe>  t ..plications  and  Conclusions 

By  means  of  Equation  (6)  a  method  of  solving  for  ci 
between  mandatory  radiosonde  levels  was  established.   The 
resulting  value  of  cj:  was  there  the  mean  value  for  the  layer 
of  thickness  A2[:2it|-  Z;  •   However  c^  turns  out  to  be 
somewhat  variable  with  height.   Thus  in  Figures  2  through 
9,  all  graphs  of  c^  vs  height  have  been  drawn  as  continuous 
functions.   This  means  that  Equation  (4)  may  be  reconsidered 
from  the  viewpoint  that  c^  be  variable  with  height  in  any 
layer.   In  this  case,  however,  the  mean  value  c±   between 
standard  levels  shall  be  the  same  as  that  computed  by 
Equation  (6). 

In  the  field  of  radar-propagation,  one  frequently  must 
deal  with  the  ray-curvature   K  =  d/3-^js  given  by 

K  =  /o"fe  (|£  +  .018)   c^  6  (9) 

where  12    is  the  ray-direction  relative  to  the  spherically 
horizontal  level  I  3  \  .   Up  to  several  years  ago,  it  was 
customary  to  use  as  a  standard  &N/$z.   a  value  characteristic 
of  the  low  troposphere,  namely  d^/dz   ~  "  H  (  °0  where  a 
is  the  radius  of  earth.   We  know  now  from  the  work  of  Bean 
and  Thayer  and  others  (including  this  work)  that  \^N/sz)d^g 
is  negative  but  decreases  in  magnitude.   Thus,  by  Equation 
(5)  with  C-L  and  N  both  decreasing  with  height  in  the  tropo- 
sphere,  ^^^ctzl   would  be  relatively  small  just  below  the 
tropopause. 

Consider  still  another  refinement  in  Equation  (4), 
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namely  that  c^  is  empirically  given  as  a  function  of  height 
of  the  type  shown  in  Figures  2  through  9,  inclusive,  for  the 
localities  considered.   Then  logarithmic  differentiation  of  Equation 


(4)  gives     /  $N 
N      dTz 


Cc  *(2  -Z.J  ^ 


&   Z  I  cH^i.(I-2i)       K,UJ 

This  expression  for  &  N/&~z_    ^as  a  correction  tacuui,  as 
compared  with  Equation  (5),  given  by  the  second  term  on  the 
right  side  of  (10).   Actually  this  term  can  be  readily  com- 
puted in  finite  difference?  from  graphs  of  Cj,  ys    log,0(Z-Z;) 
(for  example,  see  Figures  2  through  9).   The  correction 
factor  may  take  on  some  importance  in  radar-propagation  in 
the  tropopause  layer,  ii>  view  of  the  large  increase  c.   j_ 
wit!  i  Lag  there. 

Finally  it  should  be  noted  that  if  an  equation  of  the 
Bean  and  Thayer  type,  Equation  (3),  is  desired.  It  may  be 
shown  that  c  is  a  weighted-mean  of  the  individual  c.,  or 

2  r,  +  (  C    =   C„  AZ  „  +    •  •   "   +    C,  A  Z  ,        (U) 

Equation  (11)  makes  it  apparent  that  the  use  of  a  mean 
c  for  the  troposphere  is  reasonable.   However,  it  suggests 
that  one  should  use  a  different  weighted-mean  for  the  lower 
stratosphere  than  for  the  troposphere.   In  Equation  (11), 
2.^^,  is  the  height  above  the  surface,  so  that  the  averaging 
must  extend  down  to  the  earth's  surface.   Equation  (11)  may 
be  readily  solved  fox  c  using  »  .     f  :•>-.  fcype  shown  in 


22 


Figures  2  through  9,  by  selecting  the  c^-values  and  weighting 
them  by  the  indicated  thicknesses  A2j,   However,  it  should  be 
remembered  that  in  most  cases  only  half  of  the  lowest  layer  t 
A  2. ,  ,  is  shown  in  these  figures.   For  the  stations  studied 
in  this  paper,  the  Ns  -values  appropriate  to  the  Bean  and 
Thayer  exponential  form,  Equation  (3),  are  listed  in  the  first 
row  of  Table  1. 
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